Weak and Classical Solutions to Multispecies Advection–Dispersion Equations in Multilayer Porous Media

The basic model motivating this work is that of contaminant transport in the Earth’s subsurface, which contains layers analytical and semi-analytical solutions corresponding advection–dispersion equations could be derived. Then, using interface relations between adjacent layers, one can streamline study to solution initial boundary value problem for a coupled parabolic system on partitioned domains. For IBVPs, we set up weak formulations prove existence uniqueness appropriate Sobolev-like spaces. A priori estimates at different levels input data smoothness were obtained. nonnegativity preservation over time discussed. We numerically demonstrate how solve reduced truncated instead original multispecies with large number layers.